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Abstract. Sea level and surface current variability in the Gulf of St Lawrence 
have been investigated primarily using seven years of TOPEX/Poseidon altimeter 
data. An orthogonal response analysis is used to derive an annual cycle from 1-s 
altimetric data along satellite ground tracks, while simultaneously removing 
aliased residual tides and dynamic signals at alias tide frequencies. An 
examination of tidal-frequency variability points to the need for a better tide 
model for detiding altimetric data in order to study shorter (than seasonal) 
period processes in the Gulf of St Lawrence. Annual sea level amplitudes and 
phase fields are constructed from the along-track analysis results using a linear 
interpolation procedure. The altimetric annual harmonic has a magnitude of 
2—Scm in amplitude and is highest in fall. The altimetric sea level results agree 
well with independent tide-gauge data at coastal stations and can be accounted 
for mainly by steric height. Geostrophic surface current anomalies derived from 
the altimetric annual sea level anomalies are then discussed in conjunction with 
numerical solutions from a regional hydrodynamic model. Interannual sea level 
change in the Gulf of St Lawrence is also investigated from both altimetry and 
tide-gauge data, indicating a prominent gulf-wide sea level drop around 1997, 
with exact timing dependent on location. The interannual sea level variability is 
thought to be primarily associated with the Labrador Current transport 
variability (via both the Strait of Belle Isle and Cabot Strait) and the south— 
north fluctuation of the Gulf Stream position (via Cabot Strait). 


1. Introduction 


The Gulf of St Lawrence (GSL) is a semi-enclosed sea (figure 1) with high 
‘stratification in the vertical. It communicates directly with the North Atlantic 
Ocean through Cabot Strait between Newfoundland and Nova Scotia, and to a 
lesser degree through the Strait of Belle Isle between Newfoundland and Labrador. 
Both Straits admit the shelf water into the GSL (Petrie et a/. 1988, Han et al. 1999) 
but only Cabot Strait allows the penetration of the deep North Atlantic water 


(Bugden 1991). A considerable amount of freshwater enters the GSL mainly from 
the St Lawrence River year round but with notable seasonal variability 
(Koutitonsky and Bugden 1991). Therefore, the GSL as a whole behaves as a 
large estuary with significant roles of the Coriolis force and baroclinic features in 
the dynamical processes. The basic estuarine circulation patterns are modified by 
bottom topography and coastline geometry. The Laurentian Channel, with depths 
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Figure 1. Map showing the study area. Open circles are the locations of coastal tide-gauge 
stations. The thin lines stand for the TOPEX/Poseidon ascending (southwest— 
northeast) and descending (northwest-southeast) ground tracks. The 100m and 
200m isobaths are also shown. AC: Anticosti Channel; AI: Anticosti Island; CT: 
Charlottetown; EC: Esquiman Channel; JCS: Jacques Cartier Strait; LE: Lower 
Escuminac; MB: Miramichi Bay; MI: Magdalen Islands; MS: Magdalen Shallows; 
PEI: Prince Edward Island; RR: Renard River; SBI: Strait of Belle Isle; SI: Sept-isle. 


of over 300m throughout its 1100 km length from the shelf edge to the St Lawrence 
Estuary, separates the southern Gulf (including the Magdalen Shallows and 
Northumberland Strait between the New Brunswick and Prince Edward Island) 
and the northern GSL (including Jacques Cartier Strait, Anticosti Channel and 
Esquiman Channel). In addition to the buoyancy forcing induced from the 
freshwater input, the dynamic processes of the GSL are also driven by 
meteorological forcing, tidal forcing and other oceanic forcing at various scales. 
The tidal variability in the GSL is predominately forced by tides in the Atlantic 
Ocean. 

Upper-layer circulation in the GSL is generally cyclonic (Koutitonsky and 


Budgen 1991, Han et al. 1999). ‘The flow through the Strait of Belle Isle is 
dominated by the inflow of the Labrador Shelf water (Petrie e¢ a/. 1988), while the 
northeastern Cabot Strait features an inflow of the Newfoundland Shelf water and 
the southwestern Strait features the outflow toward the Scotian Shelf (Han er al. 


1999). The northeastern side of the deep Laurentian Channel 1s characterized by a 
northwestward inflow originating at the offshore end of Laurentian Channel and 
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varying at interannual scales, with the mixture of the Labrador slope water from 
the north and the warm slope water offshore. 

To date, there is little information on the application of satellite altimetry to 
oceanic circulation in the GSL. This situation is caused by particular difficulties in 
extracting oceanic signals for the GSL from altimeter data. The first difficulty lies in 
the small signal-to-noise variance ratio. The second is the complex geometry of the 
Gulf and the proximity to the land. The third problem is related to the tidal 
correction; the tidal model used for detiding altimetric data was not accurate 
enough. Besides, the ice coverage in winter results in gaps in altimetric sea level 
data. However, the on-going TOPEX/Poseidon (T/P) mission has reached an 
unprecedented accuracy of 2cm for measuring temporal variability of sea surface 
height. Most of the major tidal constituents are separable from the annual cycle and 
longer-term variability in the T/P measurements. Therefore we are encouraged to 
examine sea level and surface current changes in the GSL from T/P data. 

In this study we investigate seasonal and interannual changes of sea level and 
explore mechanisms underlying their variability in the GSL. Annual harmonics of 


sea level variation and residual tides are extracted from altimetric sea surface 
heights using a modified response analysis method. Seasonal surface current 
anomalies are then calculated under the geostrophic approximation from the 
annual sea levels. The altimetric sea level and currents are evaluated and discussed 
in conjunction with in situ data and model results. 


2. Methodology 
2.1. Altimeter data processing and analysis techniques 

We have used l-second corrected collinear T/P sea surface heights with an 
along-track resolution of ~6km produced by the NASA Pathfinder Project. 
Corrections were made to account for various atmospheric and oceanographic 
effects, such as ocean tides. In particular, a local tide model (Lambert et al. 1991) 
given on a 0).25° grid, essentially the same as Godin (1980), was used for the GSL. 
The data on three ascending and five descending tracks (see figure 1) from October 
1992 to September 1999 are chosen for further processing. The TOPEX/Poseidon 
satellite repeats its ground track every 9.9156 days. Ideally there are 276 cycles of 
data for a given track location but the actual number of data is substantially lower 
due to ice coverage in winter. A location is included for analysis only if there are 
more than 200 cycles of observations available. 

To date the marine geoid still has considerable errors at spatial scales of interest 
here. To bypass the problem we have calculated a mean sea-surface height at each 
track location from all available data, and then removed the mean from them, 
producing sea level anomalies. As a result information about the mean sea surface 
topography and associated geostrophic currents is lost. This loss of information 1s 
not a concern in this study since we are interested in annual and interannual 
variability only. 

The sea surface height anomalies are then subject to a modified response 
analysis (e.g. Cartwright and Ray 1990, Han et al. 2002) to separate an annual cycle 
from variations at alias frequencies of oceanic tides. We apply the method to the 
initially corrected T/P sea surface height anomalies to extract the annual cycle and 
to remove the remaining variability of major semidiurnal and diurnal tides. 
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2.2. Linear interpolation method 

T/P altimetry provides sea surface height measurements along ascending and 
descending satellite ground tracks (figure 1). The cosine and sine coefficients of the 
annual cycle on the track locations are spatially mapped onto a regular grid of 
0.25° in longitude by 0.25° in latitude, by an objective analysis method. The de- 
correlation scale used is 3° in longitude by 2° in latitude, in consideration of the 
cross-track spacing. We then construct amplitude and phase fields of the annual 
cycle from the interpolated cosine and sine coefficients. The same procedure is used 
to generate altimetric estimates at coastal tide-gauge locations for comparison with 
in situ observations. 


2.3. Calculation of surface geostrophic current anomalies 

Sea level anomalies on the 0.25° by 0.25° grid on 1 February, 1 May, | August 
and 1 November are derived from the altimetric annual cycle. Under the 
geostrophic approximation the sea level anomalies are used to calculate surface 
current anomalies on the four dates, nominally representing winter, spring, summer 
and fall, respectively. 


2.4. Tide-gauge data 

We obtained hourly sea level data at tide-gauge stations (see figure1 for 
locations) along the coastal boundary of the GSL from the Marine Environmental 
Data Service, the Fisheries and Oceans Canada. Monthly-mean sea levels are 
calculated from the hourly data and adjusted for local inverse barometric effects for 
the study period. Monthly mean sea level anomalies are then computed by 
removing time means from the adjusted sea levels for each station. A least squares 
regression 1s performed to extract annual and semi-annual cycles from the monthly 
anomalies. The residual monthly data with the annual and semi-annual cycles 
removed are then used to generate seasonal sea level anomalies that are further 
smoothed with a temporal 5-point moving filter. 


3. Variability at the /, and K, tidal alias periods 

With the repeat cycle of approximately 10 days for the T/P altimeter, semi- 
diurnal and diurnal tides in the T/P data have much longer alias periods (Han et al. 
1996), for example, 62 days for M> and 173 days for K,. Therefore, the response 
analysis may reduce both residual tidal variability and non-tidal oceanic features at 
the alias periods. Fortunately, the alias periods for the major tides are sufficiently 
shorter than half a year (except for the alias period of K,) for the response analysis 
to extract the annual and interannual variability of interest in this study, given the 
data duration of longer than seven years. 

As introduced in §2.1, the tide model of Lambert et al. (1991) for the GSL is 
used to remove oceanic tides from the altimeter data in this study. It is expected 
that a local model is more accurate than a global tide model for this semi-enclosed 
sea. Here we are interested in how good this local model really is. Therefore, this 
section examines the variability at the M, and K, alias periods, which is removed 
from the tidally corrected altimeter data by the modified response analysis. 

The semi-diurnal tide system features amphidromes near the Magdalen Islands 
and in the western entrance of the Northumberland Strait (Koutitonsky and 
Bugden 1991), with amplitude increasing from about 30cm near Cabot Strait to 
130cm in the western Gulf. Figure2 shows the spatial distribution of residual 
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Figure 2. Amplitudes (cm) of variations at the TOPEX/Poseidon alias periods of (a) M> and 
(6) K, derived from the TOPEX/Poseidon altimetric data in the Gulf of St Lawrence. 


amplitudes estimated at the M, and K, alias periods from the corrected altimeter 
data. At the alias period of M>, (figure 2(a)), the amplitude of the fluctuation 
increases from 3cm at Cabot Strait and the Strait of Belle Isle to about 8cm near 
Miramichi Bay and about 10cm west of Anticosti Island. The large residual tide 
suggests an alternative model with higher accuracy should be considered in detiding 
altimetric data in the GSL. 

The diurnal tidal height variability is generally weaker than the semi-diurnal 
tides in the GSL, ranging from ~10cm on the east to ~20cm on the west and 
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Figure 3. Annual harmonics of sea level from (a) TOPEX/Poseidon data, (5) steric height 
and (c) amplitude and phase differences (steric height minus TOPEX/Poseidon) in the 
Gulf of St Lawrence. The amplitude contour intervals are 0.5cm. The phase indicates 
the time (year day) when the annual sea level is highest. 


south. In the tidally corrected T/P data, the response analysis indicates residual 
variability at the K, alias frequency is relatively uniform in space, at about 2cm 
(figure 2(b)). Any substantial variability at the alias period of K, can be due to 
either the tide correction error or the semi-annual variability. The derived 
variability at the K, alias period in the GSL seems more likely related to the tidal 
fluctuation, since an analysis of the tidal gauge data indicates smaller amplitude for 
the semi-annual cycle. 
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Figure 3. (Continued) 


Although the tide model of Lambert et al. (1991) is believed to be better than 
global tide models for the Gulf, the present analyses indicate its substantial errors 
for M, (figure 2(a)). The tidal errors are not a problem for the present investigation 
of the annual and interannual variability, but can be for a study of dynamical 
processes at semi-annual or shorter scales. 


4. Seasonal sea level and surface current variability 
The response analysis extracted the annual cycle of sea level while 
simultaneously removing sea level variability at the tidal alias frequencies (discussed 
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Figure 3. (Continued) 


in the preceding section) in the T/P sea level time series. In this section, we examine 
amplitude and phase of the annual cycle of altimetric sea levels in the GSL, 
compare the altimetric results with coastal tide-gauge measurements, and explore 
altimetric sea surface current variability. 

The magnitude of the annual cycle varies from 2—5 cm, smaller at the confluence 
region of the Laurentian Channel, Esquiman Channel and Anticosti Channel 
(figure 3(a)). The phase (indicating the time of annual maximum sea level) pattern 
shows substantial changes with a general lag (about two months) from the western 
Gulf towards the eastern Gulf. Overall the highest sea level associated with the 
annual harmonic occurs in September—November. 
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The amplitude of the annual cycle is ~5cm in the Magdalen Shallows and 
Northumberland Strait (figure 3(a)). The sea level is highest in late summer/early 
fall in the west and in fall in the east. This seasonal variation is probably associated 
with steric effects due to thermal expansion in response to the seasonal surface 
heating and due to the seasonal evolution of seaward freshwater advection from 
late spring through late summer. 

A sea level minimum is located in the Laurentian Channel near Cabot Strait. 
Another sea level minimum is situated over the confluence zone of the Laurentian 
Channel, Esquiman Channel and Anticosti Channel. The former is qualitatively 
consistent with the seasonally persistent circulation pattern of an inflow into the 
Gulf of St Lawrence on the Newfoundland side of the Laurentian Channel and an 
outflow from the GSL on the Nova Scotia side (Han et al. 1999). The latter may be 
associated with a topographically controlled eddy (with little seasonal change in 
strength) over the confluence zone. 

The sea level variability increases northeastward toward the Strait of Belle Isle. 
The phase shows a sea level peak in late fall. Earlier tide-gauge data showed 
significant seasonal sea level differences across the Strait of Belle Isle, with a 
maximum in fall and minimum in summer. Estimates from water mass analysis, 
moored measurements (Petrie et a/. 1988) and numerical modelling (Han et al. 
1999) indicate stronger inflows in the fall/winter and weaker inflows in summer. 
Hence, the sea level variability of the northeastern GSL appears closely related to 
inflow through the Strait of Belle Isle from the Labrador Shelf. 

The sea level in Jacques Cartier Strait shows a substantial increase westward, 
which is probably associated with amplification due to flow convergence in the 
Strait. The phase pattern (the high sea level in Jacques Cartier Strait in early fall 
versus that in the Strait of Belle Isle in late fall) indicates that the effect of the local 
heating/cooling dominates over that of the advective inflow from the Strait of Belle 
Isle. 

Comparisons of the T/P derived annual cycle with tide-gauge measurements 
(table 1) show overall good agreement. The root mean square (rms) amplitude 
and phase differences are 0.4cm and 11 days, respectively. The good agreement 
demonstrates the ability of the T/P altimetry’s observing seasonal sea level 
variability in the GSL. It also suggests that in spite of data unavailability in the ice 
season the modified response analysis is able to reliably extract annual variability, 
taking advantage of the interannual variation of the ice season and the long data 
span. The remaining discrepancy may be attributed to altimetric data degradation, 
spatial sparsity and extrapolation. 

We calculated steric heights from the climatological seasonal-mean density fields 
of Han et al. (1999) and derived annual harmonics of steric height (figure 3(b)). In 


Table 1. Annual harmonics at selected tide-gauge stations. Amplitudes are in cm and phases 
indicate the year day (month) when the annual sea level is highest. 


Tide gauge TOPEX/Poseidon 
Amplitude Phase Amplitude Phase 
Sept-iles 4 253 (September) 3.8 248 (September) 
Renard River 5.) 257 (September) 3.6 255 (September) 
Charlottetown 4.3 317 (November) 5 309 (November) 
Lower Escuminac 3.5 268 (September) 27 287 (October) 
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general, the altimetric annual cycle agrees fairly with the steric height (figure 3(c)). 
The mean amplitude of the T/P annual cycle is 3.3cm, with a mean difference of 
0.8cm (steric height minus altimetry) and an rms difference of 1.3cm. The mean 
and rms phase differences are —11 (steric height minus altimetry) and 25 days, 
respectively. The fair agreement supports the preceding discussions that the annual 
sea level variability as observed by T/P altimeter can be largely accounted for by 
steric height variability, which is mainly associated with local heating/cooling, river 
run-off, and precipitation/evaporation. 

The geostrophic surface current anomaly on | February, | May, 1 August and 
1 October is derived from the annual cycle of the altimetric sea surface height. The 
T/P data do not provide accurate information on the mean current at the spatial 
scales of interest due to the uncertainty associated with the present marine geoid. In 
the future, a sufficiently improved geoid from the Gravity Recovery and Climate 
Experiment mission can be used. As an alternative, we use an independent 
climatological-mean surface current field from a diagnostic finite element model for 
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Figure4. (a) Climatological from Han et al.’s (1999) 
diagnostic model in winter, spring, summer and fall. (6) Annual variation of sea- 
surface currents, constructed from Han ef al.’s (1999) annual means and the 
altimetric anomalies on 1 February, 1 May, 1 August and 1 November. 
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the GSL (Han et al. 1999). The model circulation is driven with climatological 
seasonal density fields, observationally based sea level across the Strait of Belle Isle, 
and steric height conditions at the other open boundaries (see Han et al. (1999) for 
details). We have interpolated the model seasonal-mean (January—February, April— 
May, July-August and October-November nominally for winter, spring, summer 
and fall, respectively) surface currents to the 0.25° by 0.25° grid (figure 4(a)) using 
the linear interpolation method (see §2.2). The model results provide detailed 
seasonal-mean surface circulation patterns on the gulf-wide scale, producing 
dominant cyclonic circulation with inflows through the Strait of Belle Isle and on 
the northeastern side of Cabot Strait, and outflow on the southwestern side of 
Cabot Strait. These inflows circulate on the northern Gulf towards the western 
Gulf, mix with St Lawrence River run-off, and move seaward on the southern Gulf 
(Han et al. 1999). The apparent flow crossing some of the northern coast is an 
artificial outcome from the sub-optimal linear interpolation method. 
Seasonal-mean current fields (figure 4(b)) were also constructed by adding the 
altimetric surface current anomalies to the model mean currents from the four 
seasonal solutions. The overall GSL circulation is strongest in fall and weakest in 
spring. The seaward flow along the western side of Laurentian Channel and on the 





(b)52° 52° 
51° 51° 
50° 50° 
49° 49° 
\\ 
Sy 
48° : nt eo 48° 
Aina A pe 
t 
47° 47° 
46° 46° 
45° 45° 
-68° -66° -64° -62° -60° -58° -56° -68° -66° -64° -62° -60° -58° -56° 
52° 590 
51° 51° 
50° 50° 
49° 49° 
\ es 

48° : > oN 48° = 4 4 \ a 7 iN 

Ke Wyse ayes a Sy 

a | SUES ESSN &* 


Wee > 
: NAGE AA 
14 = yj 
46° ee ¥ *~- 


45° 
-68° -66° -64° -62° -—60° -—58° -56° -68° -66° -64° -62° —60° -58° -56° 





Figure4. (Continued) 
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southern Gulf appears to be larger in summer/fall. The seasonal current variability 
derived from the T/P altimetry is generally consistent with Han et al.’s (1999) model 
results (figure 4()) in many aspects, for example, a stronger seaward flow over the 
southern Gulf in summer and a larger outflow on the western Cabot Strait in fall. 
However, there are notable discrepancies in the detailed circulation features 
(figure 5). Quantitatively, the rms current values in the GSL are calculated to be 7.9, 
7.7, 7.9 and 8.0cms ! for winter, spring, summer and fall respectively; while the 
corresponding model-altimetric rms differences are 3.8, 3.8, 4.1 and 4.4cm a 
respectively. Note that the altimetry results provide annual harmonics of 
geostrophic current estimates for the period from 1992 to 2000, and the model 
solutions represent climatological bimonthly realization diagnosed from density 
data. Sparsity and low signal-to-noise variance ratio of altimetric data and 
inadequency of the hydrodynamic model are also factors that may have contributed 
to these discrepancies. 
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Figure5. Figure 4(a) minus figure 4(b) for the respective seasons. 
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5. Interannual sea level variability 

Interannual sea level variability is investigated by examining seasonal-mean sea 
level anomalies, after the annual cycle and tidal residuals (including the semi-annual 
cycle; see §3) are removed from altimetric data. The seasonal-mean sea level 
anomalies are calculated and smoothed with a temporal five-point moving filter for 
three regions: the northern Gulf (figure 6), the Laurentian Channel (figure 7), and 
the southern Gulf (figure 8). 

The altimetric sea level in the GSL was generally higher during the period of 
1996-1997, depending on locations (figures 6-8). The interannual range amounted 
to 5-10cm. There were notable regional differences in terms of magnitude and 
timing. The sea level drop occurred first on Track 113 in the northeastern Gulf at 
the beginning of 1996, proceeded onto Track 075 in summer and Track 033 in fall. 
In the Laurentian Channel, the sea level reached a maximum first at Track 037 in 
spring 1996, at Track 122 in summer, at Track 033 and 126 in fall. In the southern 
Gulf, the sea level peaked at the end of 1996. 

Pressure-adjusted tide-gauge data at three stations in the southern Gulf for the 
same period (figure 9, dashed lines) exhibited similar interannual variations but with 
an additional rising trend. The T/P altimeter measures a geocentric height, while the 
tide gauge measures the sea level relative to the local seabed. An early study 
(B. Petrie, personal communication, 2000) found that the sea level increase as 
observed by tide gauges along the northeastern American coast in the past century 
could, to a significant extent, be attributed to the post-glacial rebound (PGR) of the 
seabed. To account for the post-glacial rebound effect, we have de-trended the tide- 
gauge data (figure9, solid lines) by linearly fitting the data and subsequently 
removing the fitted line. We further made a quantitative evaluation by comparing 
the seasonal-mean sea level anomalies averaged for the three tide-gauges and the 
T/P data on Track 033 and 126 in the southern Gulf (figure 10). We can see from 
figure8 that the de-trended tide-gauge data are in better agreement with the 
altimetric observations both in magnitude and in phase. The correlation coefficient 
is 0.76 and the rms difference 2.0cm but without de-trending, the corresponding 
values are 0.58 and 2.6cm. The trend of sea level changes observed by the tide 
gauges is estimated to be a rise of 0.8cmyear |, much larger than 0.2-0.4cm year ! 
from the PGR model (Tushinghan and Peltier 1991) and from earlier analyses of 
longer tide-gauge data (e.g. Douglas 1991). The tide-gauge results are sensitive to 
the length of record considered (B. Petrie, personal communication, 2000). 

The interannual sea level variability in the GSL can be caused by a number of 
factors, of both oceanographic and atmospheric origins. The GSL circulation is not 
only an estuarine response to the freshwater run-off from the St Lawrence River 
but also an integrated part of the northwestern North Atlantic coastal currents. The 
Labrador Shelf water enters the northeastern Gulf through the Strait of Belle Isle 
(Petrie et al. 1988) and has a strong seasonal cycle. The Labrador Slope water 
enters the Gulf through Laurentian Channel. Changes of water properties and 
volume transports in the Labrador Current can have effects in the GSL at both 
seasonal and interannual scales. The north-south movement of the Gulf Stream 
position can affect coastal sea levels. It has been found that a northward movement 
of the Gulf Stream corresponds to a lower coastal sea level over the Scotian Shelf 
and in the Gulf of Maine (Han 2002). Besides air temperature, wind and 
precipitation are factors impacting heat and water fluxes across the sea surface 
within the GSL. They affect water properties and steric heights directly and 
dynamic sea level and currents indirectly. 
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Figure 6. Interannual sea level variability from TOPEX/Poseidon data in the northeastern 
Gulf. The thick track segments in figure 1 distinguish locations over the Laurentian 
Channel from those over the northern Gulf and over the southern Gulf. (a) Track 
113, (b) track 075, (c) track 033, (d) track 037, (e) track 122. 


The consequence of the air temperature variability is assessed using air 
temperature records at Magdalen Island. High air temperature increases the vertical 
length of the upper water column through thermal expansion. We calculated 
seasonal-mean temperature anomalies from monthly air temperature data, applied 
a five-point moving filter to the seasonal anomalies, and finally generated annual 
anomalies (figure 11(a)) from the filtered seasonal anomalies for the period from 
1993 to 1999. The air temperature data indicated substantial increases in the early 
and late 1990s, becoming relatively steady in 1996-1997 (figure 11(a)). It follows 
that the altimetric and tide-gauge sea level decline, which started in late 1996 and 
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Figure 7. Same as figure6 but in the Laurentian Channel. (a) Track 037, (b) track 122, 
(c) track 033, (d) track 126. 
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Figure 8. Same as figure6 but in the southern Gulf. (a) Track 033, (b) track 126. 
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Figure9. Seasonal-mean sea-level anomalies from coastal tide-gauge data at (a) Renard 
River, (b) Lower Escuminac and (c) Charlottetown before (dashed lines) and after 
(solid lines) de-trending. 


early 1997, appears not to be a consequence of the air temperature variation over 
the GSL. 

The effect of the St Lawrence River run-off is evaluated by examining its low 
frequency variability. The run-off can significantly impact the salinity distribution 
in the GSL, especially in the southern Gulf and in the upper water column (Han 
et al. 1999). High run-off reduces salinity and hence is expected to increase steric 
height. We computed seasonal-mean run-off anomalies from Bourgault and 
Koutitonsky’s (1999) monthly St Lawrence River run-off at Quebec City, then 
applied a five-point moving filter to the seasonal anomalies, and finally produced 
annual anomalies (figure 11(5)) from the filtered seasonal anomalies for the period 
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Figure 10. Seasonal-mean sea-level anomalies from coastal tide-gauge data (dashed thick 
line: before de-trending; solid thick line: after de-trending) and altimetry (thin line) 
averaged in the southern Gulf of St Lawrence. 
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Figure 11. Annual anomalies of (a) air temperature at Magdalen Islands and (5) run-off of 
the St Lawrence River at Quebec City. 


from 1993-1999. As can be seen from figure 11(5), the St Lawrence River run-off 
had substantial interannual variability in the 1990s, highest in 1997 and lowest in 
1999. Nevertheless the sea level in the Gulf of St Lawrence started to fall in late 


1996 and early 1997. ‘Therefore the St Lawrence River run-off appears unlikely to 








We conjecture that the interannual sea level variability is primarily related to the 
variation of the Labrador Current transport and the fluctuation of the Gulf Stream 
position. The sea level change was initially associated with the change of the 
position of the Gulf Stream. In 1993-1994 when the North Atlantic Oscillation was 
very strong and the Gulf Stream was in its most northern position (figure 12(a)), the 
sea level in the GSL, similar to that over the Scotian Shelf (Han 2002), was 


relatively low. When the Gulf Stream retreated southward in 1995-1996 the sea 
level in the GSL started to rise. The trend continued until a pulse of the colder 
Labrador Current water intruded into the GSL. At first we saw the effects of the 
colder water intrusion in the northeastern Gulf in 1996, possibly due to a stronger 


inflow through the Strait of Belle Isle as a result of an increased Labrador Current 
transport (figure 12()). The gulf-wide influence did not occur until the Labrador 


Current water pulse entered the Gulf through the deep Laurentian Channel. The 
process is similar to the sea level fall in 1997/1998 on the Scotian Shelf and in the 
Gulf of Maine (Han 2002). They reported that the influence of the cold Labrador 
Current water pulse occurred in the outer Laurentian Channel between the 
Newfoundland and Scotian Shelves in late 1996. The present study revealed that 
the gulf-wide sea level drop happened in early 1997. Therefore the variation of the 
mixed slope water at the shelf break propagated along the Laurentian Channel into 
the Gulf, as stated in Bugden (1991). Afterwards, though effects of the Labrador 
Current water pulse faded out, the onshore movement of the Gulf Stream kept the 
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Figure 12. (a) Annual anomalies of the Gulf Stream northern wall position (K. Drinkwater 
and R. Pettipas, personal communication, 2001). (6) Spring/summer baroclinic 
transport anomalies of the Labrador Current (Han and Tang 1999). Positive 
anomalies indicate larger southward transport. 


GSL sea level remaining low. Therefore, the interannual sea level variability 
appears to be controlled by the interplay between the baroclinic Labrador Current 
transport and the proximity of the Gulf Stream to the eastern Scotian and 
southwestern Newfoundland Shelves. 


6. Summary 

We have investigated annual and interannual sea level variations in the GSL, 
using the T/P altimeter observations in conjunction with tide-gauge data in the 
1990s. Altimetric geostrophic current anomalies are also calculated and described 
together with the mean currents from a diagnostic finite-element model (Han et al. 
1999). Although the tide model (Lambert et al. 1991) is not perfect, its residuals, 
which are aliased by the T/P sampling scheme and can be removed by the present 
response analysis method, seem not critical to the investigation of the annual and 
interannual changes. 

The present analysis reveals a substantially variable spatial distribution of the 
annual sea level cycle, with the amplitude changing from nearly 2cm in the central 
GSL to about 5cm in the southern GSL and the sea level peaking in the 
northwestern GSL in late summer and in the eastern Gulf in early fall. The phase 
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The good agreement between the altimetric results and tide-gauge observations 
suggests that the modified response analysis is capable of accurately recovering 
annual variability, overcoming data gaps in the winter ice season and the 
unfavourable signal-to-noise ratio. 

The seasonal variation of geostrophic surface currents is derived from the 


annual cycle of altimetric sea levels. The T/P results indicate that an overall cyclonic 
GSL gyre is strongest in fall and weakest in spring. The altimetry-derived current 
variability is generally and roughly consistent with the seasonal change of the model 
GSL surface circulation (e.g. Han et al. 1999). 

The striking feature of the interannual variability in the GSL for the study 
period is the sea level rise prior to 1996, and fall in 1996 and 1997. The magnitude 
of the sea level fall 1s at S—-10cm, larger than that of the annual cycle. The altimetric 


variations are consistent with tide-gauge observations. This interannual change 
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